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Abstract 



We report on a study of interaction effects in the tunneling density-of- 

states of a disordered two-dimensional electron gas in the strong magnetic 

field limit where only the lowest Landau level is occupied. Interactions in the 

presence of disorder are accounted for by performing finite-size self-consistent 

Hartree-Fock calculations. We find evidence for the formation of a pseudo-gap 

with a tunneling density-of-states which vanishes at the Fermi energy. 
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In the absence of a magnetic field, interaction effects in disordered electronic systems 
cause dramatic reductions in the tunneling density-of-states (DOS), gri^), at energies near 
the chemical potential fi. In two-dimensions in the weak disorder limit the relative suppres- 
sion of QrifJ') increases logarithmically with decreasing temperature At sufficiently low 
temperatures, where the weak-disorder approximation fails, all states are localized. For lo- 
calization lengths small compared to the distance between electrons (the classical limit) the 
Coulomb gap theory of Efros and Shklovskii ^ is expected to apply. In the strong mag- 
netic field limit, the one-body electronic spectrum in the absence of interactions and disorder 
consists of macroscopically degenerate Landau levels and the behavior of gri^) promises to 
be even richer. Weak-disorder theory gives the same temperature dependence as at 
-6 = 0, although because of the absence of a small parameter it cannot be formally justified 
even at high temperatures. At very low temperatures all states in disordered samples are 
expected to be localized except at a single energy within each Landau level and Polyakov 
and Shklovskii ||^ have recently argued that the low-temperature transport properties 
of disordered samples in strong magnetic fields provide indirect evidence for Coulomb gap 
formation. Direct experimental evidence of strong suppression of gri^) for e near /i comes 
from the recent experiments of Ashoori et al. and Eisenstein |Q] et al. In this Letter we 
report on the first study of the suppression of gri^) in a disordered two dimensional gas of 
quantum electrons at strong magnetic fields. 

We restrict our attention here to the strong magnetic field limit where the Landau level 
filling factor, u = N/N^ < 1, the electrons are completely spin-polarized, and Landau level 
mixing by interactions or disorder can be neglected. (Here = BL^Ly/ $o = L^Ly/ (27r£^) 
is the number of states per Landau level, $0 = hc/e is the magnetic flux quantum and 
N the number of electrons.) Coulomb interactions are treated self-consistently within the 
Hartree-Fock approximation (HFA) . As we discuss further below the HFA becomes exact in 
the classical limit of strongly localized states and is generally expected to be accurate when 
the disorder energy scale is larger than the interaction energy scale. Calculations have been 
completed for u = 1/5 and for z/ = 1/2. In the former case the states near the Fermi level 
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in the absence of interactions have a locahzation length comparable to the distance between 
electrons and we find clear evidence of the formation of a pseudogap where the DOS at the 
Fermi level vanishes in the thermodynamic limit. As far as we are aware ours are the first 
calculations which demonstrate the occurrence of a Coulomb gap away from the classical 
limit. We should expect to find qualitatively similar results at other values of z/, except 
possibly at z/ = 1/2 (see below). 

We perform our calculations in the Landau gauge (A = (0, Bx, 0)) and apply quasiperi- 
odic boundary conditions to the HF single-particle orbitals inside a rectangle with sides 
Lx and Ly. The basis states used to represent the HF Hamiltonian are related to elliptic 
theta functions and can be labeled by a set of guiding centers \Xi), i = 1, ■ ■ ■ , inside the 



fundamental cell of the finite-size system [10|. The finite-size system HF calculations in the 



absence of disorder have been described in detail previously |jTT|. An important simplifica- 



tion related to the strong magnetic field limit [|T2[ is the fact that both Hartree and Fock 



potentials can be expressed as a functional of the electron density: 

'iqx{Xi + X 



{X^\VHF\X,) = J2Aiq}exp 



6{j,i + my)UHF{q) (1) 
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Here the sum over q is over the discrete set of wavevectors consistent with the boundary 
conditions and 6{n',n) is 1 if n' = n{modN^) and otherwise. The quantity A(g) is 
proportional [|lT| to the Fourier component of the charge density at wavevector q and is 
calculated from the eigenvectors of the HF Hamiltonian by 

m = ]f E E +m,)exp [^M^^^±^j ^(x,,|a)(a|X,). (2) 



a=l 



In Eq. (0) the sum over a is over the N lowest energy eigenvectors of the HF Hamiltonian. 



Unpiq) includes both Coulomb and exchange interactions and needs to be evaluated |TT 
only once for each finite system size. UHpiq) is proportional to j kI {k is the dielectric 
constant) which is the characteristic interaction energy scale for interacting electrons in the 
lowest Landau level. 

We have used a model disorder potential consisting of Nj 5-function scatterers with 



strength uniformly distributed between —A and A and scattering centers uniformly dis- 
tributed inside the fundamental cell of the finite-size system. The disorder-averaged DOS 



in the absence of interactions is known exactly for this model ||T3[ and similar models have 
been used previously to study localization of non-interacting electrons in strong magnetic 
fields [0]. For a given disorder realization the eigenstates \a) and eigenvalues of the HFA 
satisfy the following equation: 

Y,{X,\HHF\X,){X,\a) =Y^\{Xj\Vj\X,) + {X^\Vhf\X,)\ {X,\a) = (3) 

i i 

For each disorder realization Eqs. @, ([^), and @ were solved by iterating to self-consistently 
for six different values 7 = {e^/ntj/V where F = [y^N ij P'LxL^^I'^ is the characteristic 
disorder potential energy scale. (For this model the self-consistent Born approximation 
DOS is non-zero for |e| < (2/37r)^/^r ~ 0.46r.) The HF equations were solved first at 
7 = where no self-consistency is required, and then at 7 = 7j = O.lz, % = l,---,5. 
For 7 = 7j the iteration is started using the self-consistent charge density obtained at 
7 = 7j_i. The average change in the charge- density Fourier components at iteration j, 
[Y^ql^j+iio) ~ ^jiQ)\]X^'^ = 5, was evaluated to test for convergence. The calculations were 
considered to be converged once 6 < 10^^. All the numerical results reported below were 
obtained for Nj = 5N^. Square {L^ = Ly = L) systems of various sizes were considered for 
both z/ = 1/5 and u = 1/2. 

For V 01 1 — V very small the classical limit is approached. In this limit the HF energy 
for the orbital at disorder site i is 

E. = c^.+ Y.mvm - n^^mm^ (4) 

where 0j and fi are the disorder-potential energy at site i and the occupation number at site 
i. If the overlap between different orbitals can be neglected the exchange matrix element in 
Eq. (|^) is zero unless i = j and we obtain Ei = (pi + J2j^iVijfj where Vij = {ij\V\ij). In 
the localized limit the only effect of exchange is to cancel the self-interaction energy of each 
orbital. This is just the classical expression for the energy of an electron at site i used in 



previous investigations of the Efros-Shklovskii Coulomb gap [^]. Thus the HF approximation 
captures all the physics of the classical limit but is able to include quantum effects without 
additional difficulty. 

Typical numerical results are shown in Fig. (0) for z/ = 1/5 and in Fig. (|^) for z/ = 1/2. 
In the HFA the tunneling DOS is given by the spectral density of HF eigenvalues. We plot 
the total number of eigenvalues which occur in a range of width AE around the energy 
E, Dae{E), for all disorder realization studied. (Note that the histogram box widths used 
to obtain these results are several times larger than the typical energy level spacing at the 
Fermi energy.) For a given system size the DOS exhibits a minimum which is tied to the 
Fermi level and which deepens and broadens as the interaction strength increases or as the 
disorder decreases. For a given ratio of interaction strength to disorder the DOS minimum 
becomes more pronounced as the system size increases. 

In the classical limit the difference between the energy of an orbital which is occupied 
and an orbital which is empty must exceed the excitonic energy associated with the 
transfer of an electron from the occupied site to the empty site. For a finite size square system 
with periodic boundary conditions the minimum possible excitonic energy (and hence the 
minimum possible level spacing) occurs when the displacement of the charge is (L, L) / \/2 
and equals ~ 2.285e^/fi;L. This energy provides a lower bound on the energy-level-spacing 
at the Fermi level in a disordered system of classical electrons. In the thermodynamic limit 
Efros and Shklovskii have estimated that for classical electrons the impurity-averaged DOS 
at the Fermi level in two- dimensions vanishes linearly with a coefficient which corresponds 
to an energy level spacing at the Fermi level of a/^ttc^/kL ~ 2.507e^ / kL. The dependence 
of our HF spectrum DOS minimum on system size is illustrated in Fig. (|^) where we plot 
the DOS per Landau Level at the Fermi energy as a function of x = {T /N^)/{e^ / nL). x is 
proportional to the ratio of the energy-level-spacing in a finite system of electrons without 
interactions to the energy-level-spacing at the Fermi level for a disordered system of classical 
electrons with Coulomb interactions. We see from Fig. (^) that to a good approximation 
the Fermi level DOS depends only on this ratio. For large x Coulomb interactions are 
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unimportant and gri^p) approaches the known result for non-interacting electrons, even for 
quite small system sizes. For small x the numerical results suggest that gri^F) is proportional 
to X with a slope corresponding to an energy-level spacing at the Fermi level of ~ O.Se^/nL 
for z/ = 1/5 and ~ O.Qe^/nL for u = 1/2. This decrease in coefficient of the energy-level- 
spacing at the Fermi level for large L compared to the classical limit cannot be explained 
by the finite spatial extent of the quantum electronic orbitals alone @] and is a non-trivial 
quantum effect. However the dependence of the energy-level spacing on system size is 
not altered by quantum effects; we believe that these numerical results convincingly establish 
that the HF DOS vanishes at the Fermi level in the thermodynamic limit [|l^. 

For non-interacting electrons in the lowest Landau level states at the Fermi level are 
localized except (for typical disorder potentials) at z/ = 1/2. It is therefore surprising that 
our numerical results at u = 1/5 and u = 1/2 are not more different. For u = 1/2 the 
conductivity is finite and the two-particle spectral response has a diffusive form at long 



distances [|r^. To lowest order in interactions we can construct the HF self-energy from 
the non-interacting system eigenstates. Following work on interaction effects in disordered 
metals this leads [0] to a DOS change from the exchange self-energy 

D[E)^ dE ^TiDJq^q^ (r|E-e^|)V2- ^> 

Eq. is a good agreement with our numerical results at v = 1/2 far enough from the 
Fermi energy where it is permissible to treat interactions in lowest order. In Eq. (|]) we have 
evaluated the exchange self-energy using the bare two-dimensional 1/r interaction, V{q) = 
27re^/g. It is the absence of screening in the HFA which makes the DOS change increase 
more strongly as energies approach ep than the \ia{E) dependence expected for disordered 
metallic systems in two-dimensions. We expect that screening effects neglected in the HFA 
would increase the static dielectric constant, k, and change our results quantitatively but 
not qualitatively for u ^ 1/2. For u = 1/2, however, the localization length diverges and k 
may also diverge |]TB|. Further work will be required to determine whether gri^) at z/ = 1/2 
is anomalous. 



Finally we mention the downward shift in the chemical potential due to interactions 
which is visible in Fig. (|1]) and in Fig. (|^). For the case of u = 1/5 the magnitude of 
this exchange self-energy shift in the chemical potential is much smaller than the value 
V\J n /2{e^ / nt) calculated neglecting spatial correlations of the occupied orbitals and used 
in many considerations of the competition between disorder and weak interactions. This 
reduction has important implications for the compressibility of disordered electronic systems 
in a strong magnetic field and for the magnitude of exchange-enhanced spin-polarizations. 

In conclusion we have performed finite system self-consistent HF calculations for a disor- 
dered system of two-dimensional quantum electrons in the strong magnetic field limit. We 
find that the HF energy-level spacing at the Fermi level varies as c^/kL, corresponding in 
the thermodynamic limit to a tunneling DOS per area which vanishes linearly as the Fermi 
energy is approached, for all filling factors. We expect that these results apply to systems 
with sufficient disorder that the fractional quantum Hall effect does not occur. This reduc- 
tion in gri^) near the Fermi energy results from static correlations in the self-consistent HF 
ground state which limit the phase space for adding an electron near the chemical potential. 
Similar reductions in gri^) should occur for any correlated electron state, although the be- 
havior of gri^) at very low energies could differ qualitatively in the strongly correlated fluid 
states which occur in weak-disorder systems. 
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R.C Ashoori, A.L. Efros, J. P. Eisenstein, and B.I. Shklovskii. This work was supported by 
the National Science Foundation under Grant No. DMR-9113911. 
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FIGURES 

FIG. 1. Tunneling density-of-states histograms obtained from self-consistent Hartree-Fock cal- 
culations at I' = 1/5. D^e{E) is the number of Hartree-Fock eigenvalues at energies between E 
and E + ^E including all disorder realizations. These results are for = 245 and No = 90 
disorder realizations. For these plots A£' = 0.028r. 

FIG. 2. Tunneling density-of-states histograms for u = 1/2, N^j) = 288, No = 60, and 
AE = 028r. 

FIG. 3. Finite system density-of-states minima versus x. At a given filling factor results at 
different values of 7 and different system sizes collapse onto a single curve. Dae{Ef)^ / AEN^Njj 
is the density of states per flux quantum per F energy interval. The horizontal arrows at the 
right hand side of the figure indicate the exact values for the large x (no interaction) , large-system 
limit. The plotting symbols are horizontal bars (7 = 0.1), crosses (7 = 0.2), diamonds (7 = 0.3), 
squares (7 = 0.4), and octagons (7 = 0.5). The solid lines are guides to the eye. For u = 1/2 the 
plotted {N^,Nd) are (128,90), (288,60), and (450,7). For u = 1/5 they are (80,150), (125,120), 
and (245,90). {x = -f-^{2TT/N^y/^.) The dashed line is the Efros-Shklovskii expression for classical 
electrons. 



10 



